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Abstract 



Unusually large spontaneous and piezoelectric fields in the III-V nitrides 
have led to the making of an entirely new class of two-dimensional electron 
gas. Fluctuation from a perfectly periodic binary structure in highly polar 
semiconductor alloys present the same physical situation as a random distri- 
bution of microscopic dipoles. The excess dipole distribution in the barrier 
layers is evaluated by a method similar to the virtual crystal approximation. 
It is shown that the mobility of electrons in the two-dimensional electron gas 
formed in highly polar heterostructures is intrinsically limited by scattering 
from such dipoles. 
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I. INTRODUCTION 



The last decade has witnessed a formulation of the microscopic quantum theory of po- 
larization based on a rigorous definition of polarization in a periodic system [Ref . 1-3] . The 
study of polarization and its effects was hmited to ferroelectrics before the advent of III-V 
nitride semiconductors. The spontaneous and piezoelectric polarization fields of this class 
of semiconductors is imusually high compared to other 111-Vs; in fact, it has given birth 
to a new frontier of semiconductor physics. The wide band gap and strong polarization 
fields have found wide application in high power, high speed electronic and optical devices 
[Ref. 4,5]. A large amount of effort has been devoted to the fabrication and study of high 
power AlGaN/GaN heterostructure high electron mobility transistors (HEMTs) [Ref. 6]. 
Of chief importance in such devices is the electron mobility in the two-dimensional electron 
gas (2DEG) formed at the heterointerface. The strong polarization of the III-V nitride het- 
erostructures makes a clear analysis of the effect of microscopic polarization on mobility of 
considerable interest to device designers and theorists ahke. 

At high temperatures (T > lOOii') polar optical phonon scattering is the largest scatter- 
ing mechanism; the effects of the other scattering mechanisms become dominant only at low 
temperatures. The set of scattering mechanisms needed to understand the low temperature 
2DEG mobility of the AlGaAs/GaAs and Si-metal oxide field effect transistor (MOSFET) 
is fairly complete. Theoretical mobility calculations show very good agreement with ex- 
perimentally observed values [Ref 7,8]. However, the 2DEGs in III-V nitride MDHs have 
a fundamentally new origin. The 2DEG in such heterostructures can be entirely polariza- 
tion induced, as opposed to by remote doping as in AlGaAs/GaAs MDHs and gate-induced 
inversion in Si-MOSFETs. In fact the spontaneous and piezoelectric polarization is large 
enough to produce 2DEGs without intentionally doping the barrier, leading to the novel 
concept of piezoelectric doping in such systems [Ref. 9,10]. The term 'modulation doped 
hetersostructures' (MDHs) is somewhat erroneous in such devices, so we prefer to call it the 
more general HEMTs. 

Previously published calculations of mobilities in AlGaN/GaN 2DEGs have restricted 
their analysis to the set of scattering mechanisms that exists for AlGaAs/GaAs 2DEGs 
[Ref. 11,12]. This existing set of well understood scattering mechanisms is insufficient 
for polarization induced 2DEGs. It needs to be expanded by inclusion of the effects of the 
different origins of polarization induced 2DEGs. This work is directed towards this cxtcntion 
for the fundamentally new 2DEGs. In section II, we theoretically derive the transport 
scattering rates for 2DEG electrons due to dipoles. In section III, the results are applied to 
study the III-V nitride heterostructure 2DEG mobilities. In Section IV, we conclude with 
a discussion of the results , evaluation of the implications of the newly identified scattering 
mechanism and suggestion of a method to reduce it's effect. 

II. THEORY 

Scattering by dipoles and their effects on electron transport in bulk semiconductor sam- 
ples has been studied, albeit not extensively owing to it's insignificance in the non-polar Si 
and relatively weakly polar GaAs material systems [Ref. 13,14]. However, the effect of dipole 
scattering on 2DEG electron transport has not been studied to the best of our knowledge. 
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We derive the scattering rate due to dipoles for a semiconductor two dimensional electron 
gas. 

We consider the 2DEG to be perfect (i.e. ,the extent along the z direction to be zero) 
for our derivation. Extention to the more physical case of a 2DEG with finite extent along 
the growth direction involves incorporation of the relevant form factors. 

Figure [1] shows the model for the system under consideration. The dipole charges are 
separated from each other by distance d^. and the center is a distance z from the plane 
containing the 2DEG. Spontaneous and piezoelectric polarization fields (Pgp and Ppz re- 
spectively) in wurtzite AlGaN/GaN is directed perpendicular to the 2DEG plane [Ref. 15]. 
We choose the direction of the dipole to be perpendicular to the 2DEG plane to refiect this. 

The unscreened Coulomb potential at the origin due the dipole is written as 

Vunsir, z) = ■ [ , ' - ' =]. (1) 

Aireoeb + (z ^ f )2 ^r^ + (z + f )2 

Here e is the electron charge, eg is the permittivity of free space, e;, is the dielectric 
constant of the subsrate semiconductor, and r is the in-plane radius vector. 

We evaluate the Fourier transform of this potential in the wavevector (g) space Vuns{Q) = 
/Kn.(r)e^^-MV to get 

Vuns{q,z)^- (2) 

zeoCb q 

where q is the x — y in-plane wavevector. 

The Fourier component of screened and unscreened potentials are related for a degenerate 
2D EG through the Thomas Fermi approximation of the Lindhard formula for the momentum 
dependent static dielectric constant [Ref. 16]. It is written as e2D{q) = 1 + where qTF 
is the Thomas-Fermi wavevector, defined B being the effective Bohr radius in the 

semiconductor containing the 2DEG. A valley degeneracy of one and spin degeneracy of two 
is implied. The relation is 

e2D{q) 2eoe6 q + qrp 

This is the final screened potential experienced by an electron in the 2DEG due to a 
single dipole at a distance z from the 2DEG plane. 

Scattering rate from a state |k > to a state |k -|- q > is now evaluated. Born approxima- 
tion holds good for evaluation of matrix elements for dilute dipole concentrations. Transport 
scattering rate by a dilute perturbing potential in the Born approximation is written as [Ref. 
17] 

:^ = ^Zy J ^2\Vs:r{Q)\'mk + q)-E{k)]^^, (4) 



where n^mp 2D impurity concentration, k is the 2D electron wavevector before 

scattering, E{k) and E{k + q) are the electron energies before and after scattering, the S 
function is a statement of the elastic nature of the scattering, and for our case, nj^p is 
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the density of dipolcs at each Al(Ga) plane, ri^£p^ic- Vscri^) the total screened potential 
experienced by the 2DEG due to the entire distribution of dipoles. 

Figure [2] illustrates the physical location of the dipoles in AlGaN/GaN HEMTs. Their 
origin is discussed in section III. Owing to the interface roughness, there are dipoles located 
at the interface too; however, their effect on the 2DEG mobility was found to be much less 
than the far denser distribution of dipoles in the barrier. We consider the 2DEG to be 
physically located at the centroid of the spatially extending quasi-2DEG for illustrating the 
role of dipoles. A Fang-Howard wavefunction approach would yield respective multiplicative 
form factors, and is a simple extention of the theory presented here. 

The screened potential due to the distribution of dipoles in the barrier is hence given by 
a Fourier- weighted sum over all dipoles [Ref. 16] 

If we assume that the dipole distrubution on each Al(Ga) plane are completely uncorre- 
lated, the cross terms arising in the sum cancel, and we are left with a sum over different 
planes. The complex exponential can then be factored out and therefore does not contribute 
to the matrix element. For a thick AlGaN barrier, this evaluates to 



e 



2 



2e-5(^o+co) sinh{^] 



where zq is the distance of the centroid of the 2DEG from the interface [Figure [2]], and Cq 
is the separation of the planes containing the dipoles in the barrier. 

For a degenerate gas as in a 2DEG, scattering takes place mainly among electrons with 
wavevectors near the Fermi wavevector kp — -\/27m7 where Us is the 2D EG sheet density of 
carriers. So the k in the integral can be replaced hy kp- The elastic nature of the scattering 
process leads to a relation q = 2kp ■ sin(|) where 6 is the angle between the wavevectors k 
and k -|- q. Using these facts, the expression for scattering rate simplifies to 



2D 



m 



\V"'Hq)\''^£^=. (7) 



Using the screened scattering potential of the distribution of dipoles developed in Equa- 
tion [6] in Equation [7], we get the final scattering rate due to dipoles. 



III. APPLICATION TO III-V NITRIDE HEMTS 

We concentrate on AlGaN/GaN heterostructures in our analysis. Recent work shows that 
the 2DEG in such heterostructures is created by the subtle interplay of spontaneous and 
piezoelectric polarization [Ref. 15,18,19]. The ah imtzo calculations of Bernardini et a/ show 
that the III-V nitrides have unusually large polarization coefficients (an order of magnitude 
larger than AlGaAs/GaAs systems) [Ref. 15]. The 2DEG in such heterostructures can be 
entirely polarization induced. In that sense, it is fundamentally different from the 2DEG in 
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AlGaAs/GaAs heterostructures, which is got by modulation doping. It is argued that the 
mobihty in such a 2DEG should be affected by the cause of its very existence. 

For the perfectly periodic III-V nitride crystal, the microscopic picture of polarization is 
a dipole in each primitive cell aligned along the (0001) axis. The dipole moment po — e ■ do 
{do is the effective charge separation) is related to the macroscopic polarization P by the 
relation P = po/f^, where Q is the volume of the primitive cell [Ref. 20]. P is the total 
polarization, which includes the spontaneous and piezoelectric components. 

P=P3p + Ppz (8) 

A perfect binary polar lattice thus has a periodically arranged array of dipoles with 
equal dipole moments. Such a periodic arrangement of similar dipoles has a characteristic 
wavevector, and hence does not contribute to the scattering matrix element. 

However, the 2DEG in AlGaN/GaN heterostructures is confined by a barrier due to the 
undoped Al^Gai^^N ternary alloy barrier. The alloy is a disordered system with Al and 
Ga atoms arranged in a random array such that the overall composition over any plane is 
constant over Al(Ga) planes. The difTerence in spontaneous and piezoelectric polarizations 
between AIN and GaN implies that we have a dipole moment of randomly fluctuating mag- 
nitude in the barrier. We adopt a method similar to the treatment of disordered alloys by 
virtual crystal approximation to treat dipoles in disordered polar semiconductor alloys. 

We first arrive at the dipole moments in a primitive cell of coherently strained AIN and 
GaN binary wurtzite crystals. Pgp and the piezoelectric constants for both the semiconduc- 
tors were calculated by Bernardini et. al. in their recent paper [Ref. 15]. The piezoelectric 
field in a binary wurtzite primitive cell coherently strained to a x ~ y lattice constant a{x) 
from it's imstrained lattice constant oq and c{x) from cq in the z direction is given by the 
relation [Ref. 18] 

p,4,)= 2.(^4^). 1^3.-^33^], (9) 

where 631 and 633 are the piezoelectric coefficients and C13 and C33 are the elastic con- 
stants of the crystal structure. The volume of the primitive hexagonal cell is 

n{x) = ^co{x)-al{x). (10) 
Thus the dipole moment in a strained binary crystal is given by 

PdipoUx) = {Psp + Ppz{x)) ■ ^{x). (11) 

This dipole moment is calculated for both semiconductors as Pdipoie,AiN{x) and 

Pdipole,GaN{x). 

We model the disordered Al^Gai^x^ barrier as a perfect crystal superposed with a 
randomly fluctuating dipole moment at each primitive cell. Such a virtual crystal has a 
dipole moment of magnitude 

Pdipoleiav) = X ■ Pdipole,AlN{x) + {I - x) ■ Pdipole,GaN{x) ■ (12) 

The deviation from the perfect virtual crystal at all Al sites is (1 — x) • Apdipoie where 
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^Pdipole — Pdipole,AlN{x) — Pdipole,GaN{x) . 



(13) 



The deviation at Ga sites is {—x) ■ /S.pdipoie- Since tliere are x Al sites and (1 — a;) Ga sites 
on average on a Al(Ga) plane, the average randomly fluctuating dipole moment at each site 
is 

5pdipole = e- dQ = 2- X ■ {I- X) ■ \/^Pdipole\- (14) 

The absolute value is used in adding the dipole contributions since the direction of the 
dipole is immaterial in the scattering matrix element, which involves the square of the dipole 
potential. The number of such dipoles present on each Al(Ga) plane is given by 

'^'dipole = ,/q ^, (15) 

where the in plane lattice constant Qq^x) is interpolated for the alloy. 
The polarization induced 2DEG sheet density is given by the difference in polarization 
at the interface as a function of alloy fraction x [Ref. 18], 

ns{,x) = \Pp,{Al^Gai.^N) + Psp{AlxGai.^N) - Psp{GaN)\. (16) 

Note that the 2DEG concentration does not involve any modulation dopants; it is formed 
entirely to satisfy the discontinuity in the polarization fields at the interface. 



IV. RESULTS AND DISCUSSION 

The mobility inhibited by dipole scattering alone lJ-a?poie = ^^d§oie/''^* evaluated for 
different alloy compositions. The results are plotted in Figure [3]. We also plot the polar- 
ization induced 2DEG sheet density in the same figure for easy comparision. In Figure [4], 
we plot the dipole scattering inhibited 2DEG sheet conductivity given by G = efi^f^^^i^Us, 
and finally, in Figure [5] , we plot the dipole scattering limited mobility as a function of the 
2DEG sheet density. 

An expected increase in mobility with the increase in the binary nature of the alloy barrier 
is seen. The intrinsic low temperature mobility limit in the ,t = 0.1 to x = 0.4 range (which is 
typical of state of the art 111-V nitride HEMTs) is in the 400, OOOcm'^/Vs to 200, OOOcm^/Vs 
range. It is well worth noticing that this is much lower than the record low temperature 
mobilities (r^ lO^cm^/Vs) of AlGaAs/GaAs modulation doped heterostructures, and an 
order of magnitude higher than the record high mobihties in AlGaN/GaN HEMTs observed 
till date (51, 700cm'^ /Vs) [Ref. 19,21]. It hints at some more severe scattering mechanism(s) 
that determine the low temperature mobility in the III-V nitrides. 

Interface roughness was initially thought to be a mobility limiting mechanism owing to 
a) very high 2DEG sheet densities, and b) trying to fit experimental data to the existing 
set of scattering mechanisms for AlGaAs/GaAs MDHs [Ref. 12,22]. Recent experiments 
with double heterostructures point towards the contrary [Ref. 23]. The III-V nitrides is a 
fundamentally new system, and it is necessary that we address issues that make the material 
system so different from the AlGaAs/GaAs systems. 
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Dislocation scattering was identified as becoming dominant at high dislocation densities, 
which result from the lattice mismatch of the epitaxial GaN layer with the present substrates 
of choice - SiC or Sapphire [Ref. 24] . The novel method of lateral epitaxial overgrowth (LEO) 
is a promising candidate for reducing the density of dislocations in the nitrides [Ref. 25]. 
The effects of dipole scattering will be the next hurdle to overcome in pushing the mobihties 
higher. Digital alloy growth is suggested as a technique to reduce the severity of dipole 
scattering. By growing either purely Al or purely Ga layers, we use periodicity to overcome 
the scattering originating from the random nature of the alloy. However, digital alloy growth 
suffers from interdiffusion of atoms in the growth process, so dipole scattering cannot be 
completely eliminated by this method. 
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FIGURES 




2DEG 



FIG. 1. The location of the dipole with respect to the 2DEG is shown. The dipole axis is taken 
to be perpendicular to the plane of the 2DEG, keeping with the direction of the polarization field 
in the AlGaN barrier of AlGaN/GaN HEMTs. The distances used in the text in the derivation of 
the scattering rate are defined. 
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FIG. 2. The distribution of the dipoles in the AlGaN barrier is shown. Every Al(Ga) plane 
has dipoles in each primitive cell. The dipole moment at Al sites is higher than that at the Ga 
sites owing to the higher spontaneous polarization and piezoelectric constants in AIN than in GaN. 
This fluctuation leads to a random distribution of dipole moments which leads to scattering of 
the electrons in the 2DEG. The 2DEG is assumed to be located entirely at the centroid of the 
quasi-2DEG distribution for simplicity. 
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FIG. 3. Mobility of electrons in the 2DEG inhibited by dipole scattering alone is plotted as a 
function of alloy composition. Sheet density of carriers also changes with alloy composition, and 
is shown in the lower half. Dipole scattering dominates at alloy compositions in the x = 0.2 — 0.5 
range. As the binary nature of the alloy increases, dipole scattering reduces, leading to higher 
mobilities. 
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FIG. 4. Two-dimensional sheet conductivity G = efi'^^gi^Us, inhibited by dipole scattering 
alone is plotted as a function of alloy composition. 
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FIG. 5. Mobility limited by dipole scattering alone is plotted as a function of sheet density of 
carriers. The mobility reaches a minimum at 2DEG sheet density of !=a 2 x 10^^ /cm?. 



13 



